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Introduction 


From the time that wind tunnel wall interference was recognized to be 
significant, researchers have been developing methods to alleviate or account for it. 
Despite the best efforts so far, it appears that no method is available which completely 
eliminates the effects due to the wind tunnel walls. This report will discuss procedures 
developed for slotted wall and adaptive wall test sections of the Langley 0.3-m TCT to 
assess and correct for the residual interference by methods consistent with the 

transonic nature of the tests. 



WIAC Concept 

The underlying concept of both procedures is depicted below. There are two 
basic elements: the wind tunnel which generates the flow in which measurements are 
made, and the computer which now solves two related flow problems. In the full 
nonlinear correction procedure at least two transonic How problems are solved on 
the computer. The first is an equivalent inviscid tunnel flow where measured 
pressures near the wall and on the model are used as boundary conditions. The result 
of this first calculation is an equivalent inviscid model defined in terms of either its 
shape or its distribution of singularities. The second problem to be solved on the 
computer is a sequence of inviscid transonic calculations in which the equivalent 
model is used as the inner boundary condition and free-air conditions are used at the 
outer boundaries. The freestream Mach number and angle of attack are perturbed 
during this sequence in order to satisfy a best-fit criterion for the calculated model 
pressures and the measured model pressures. The two results obtained from these 
computer calculations are: corrections to the freestream conditions M and a, and 
a measure of residual interference. 
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Transonic WIAC Codes 


This WIAC concept was conceived by Kemp (ref. 1), developed into the 
TWINTAN code for slotted wall test sections (ref. 2) and extended by Kemp and 
Adcock (ref. 3) to include the effects of the tunnel sidewall boundary layer. The 
resulting code, TWINTN4, was enhanced by Green (ref. 4) to allow data on shaped 
walls to be used as the outer boundary condition for the equivalent tunnel flow 

calculation. 


• Transonic flows 

• Broad range of lift coefficient 

• Nonlinear TSDE 

• Uses measured wind tunnel data in BC s 

— Airfoil C p , C| and C d 
— Top & bottom wall C p ^ 

— Tunnel empty SWBL 8 and H 

• Three SWBL approximations 

— 2-wall (top & bottom only) 

— Barnwell-Sewall SWBL approximation 
— Murthy SWBL approximation 

• Two codes 

— Kemp’s TWINTN4 for slotted wall 
— TWNTN4A for adapted wall 



Nonlinear TSDE for WIAC 


The TWINTN4 code (and Green’s derivative) performs the flow calculations 
using the nonlinear Transonic Small Disturbance Equations (TSDE). The three 
dimensional effects of the sidewall boundary layer are incorporated into the two 
dimensional TSDE after Barnwell and Sewall (ref. 5) by the term, S. The effect of 
model aspect ratio was determined by Murthy (ref. 6) as a simple modification to the 
Bamwell-Sewall method. 


• Solves 2-D Transonic Small Disturbance Equation 
(TSDE) 
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• Three VLOR solutions 

— In -tunnel -> effective inviscid body 
— Free-air -> M cor and 
— Free-air -> interference field 
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Cartesian Grid for WIAC 


The flow field is discretized onto a Cartesian grid which is similar for both the 
tunnel flow calculation and the free air calculations. The top and bottom wall data is 
applied on grid lines included in the free air grid at the mean location of the walls. 
The data on the airfoil surface is applied at the slit on the tunnel centerline (or mean 
location of the model). The boundary condition at the inflow plane of the wind tunnel 
was left undetermined from wind tunnel data. This remaining boundary condition is 
assumed during the first pass through the correction code and approximated by 
iteration based on the difference between the computed inclination of the equivalent 
inviscid model and the geometric model according to the method devised by Gumbert 
et al (ref. 7). The first approximate iterated value is used in the second pass through 
the correction code; a third pass may be required. 
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WIAC Procedure 


In order to more easily apply the individual codes to the data thev were 
tncorporated into a procedure by Gumbert, et al (refs. 7 and 8) to pass data from one 
code to the next in a somewhat automated manner. This procedure was hist used for 

making correcnons to several data sets in order to validate the procedure and the 
individual codes. 
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Validation of WIAC Procedure 


The validation of the WIAC procedures (refs. 4 and 9) was accomplished by 
two types of comparisons. First, the corrected data was compared to the best available 
independent free-air computer code solutions. For the earlier slotted wall data 
comparisons (ref. 9), solutions from the conservative, transonic, full-potential equation 
(with viscous/inviscid interaction) GRUMFOIL code (ref. 10) were used. For the latter 
adaptive wall data comparisons, solutions from a Navier-Stokes code (ref. 11) were 
used. Second, the corrected data from several tests of the same airfoil shape were 
compared for consistency. 


• Comparison of Corrected Data With Indapandant 
Free-Air Calculations 

• Consistency of Corrected Data From Separate 
Tests of Common Airfoil 


Test Data Matrix 


The three tests of the CAST 10 airfoil under consideration here were conducted 
in the 8- X 24-inch slotted wall test section over a span of several years. Two of the 
tests were conducted using a six-inch-chord model. During the period between the 
two tests, several changes were made to the test section toaccommodate different 
instrumentation and flow visualization techniques. The other test used a three-inch- 
chord model. It was the only non-six-inch chord model tested in the 8- X 24-inch 
slotted wall test section. More specific information about the tests can be found in 
references 12 through 15. The figure shows the ranges of Mach number and Reynolds 
number over which the three tests were run. The WIAC procedure was applied to data 
for the three tests at those conditions which are similar for all three tests. These 
eleven common points are denoted as 9 in the figure below. 
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Preprocessor Plots of Airfoil C p 

The first step in the WIAC procedures is the preprocessor code where the 
primary function is to select only the pertinent information from the data tapes and 
generate an input file for TWINTN4. In the process it generates plots of the 
uncorrected data which are to be used as inner boundary conditions for the WIAC 
code solutions. Shown in the figure are the uncorrected pressure coefficient 
distributions on the model for each test at nearly the same conditions: Re c = 15 million, 
Mref = 0.765, and C L = .55. 


M re , = 0.765, C L = 0.55, Re c = 15x10 s 
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Preprocessor Plots of Wall C p 

The preprocessor also generates plots of the pressure coefficient distribution on 
the center slats of the top and bottom walls as shown in the figure. The Circled points 
in the bottom figure indicate data that was conspicuously inconsistent. The data point 
over the leading edge was removed and the data point ahead of the model was 
modified as shown by the filled square symbol. These are the data to be used as outer 
boundary conditions for the WIAC code solutions. 


M ref = 0.765, C L = 0.55, Re c = 15x10® 
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Uncorrected Lift Curves 


The correction to the angle of attack can best be shown in plots of the lift 
coefficient versus angle of attack. Shown in the figure is the comparison of the 
uncorrected lift curves for the three tests at M ref =0.73 and Re c = 10 million. For 
comparison, the results from GRUMFOIL are shown. The data from the three tests are 
quite scattered and each shows a different slope. 


M ref = 0.730, Re c = 10x10 6 - 
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First Pass WIAC Lift Curves 


The results from the first pass through the correction code, TWINTN4, are 
shown below. The lift curve slopes seem more consistent between the three tests, yet 
there is an unresolved shift between the data sets and with respect to the GRUMFOIL 
curve. 


M ref = 0.730, Re c = 10x10 6 
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Second Pass WIAC Lift Curves 


The results from the second pass through the correction code are shown below. 
All three tests show good agreement over the low lift range and the comparison with 
the independent free-air code is good. However, the data from test 169 tend to be 
inconsistent at moderate lift and all three data sets show different behavior near 
maximum lift. The early breakdown of the test 169 data and its correction may be due 
to the known inaccuracy of the top wall pressure data in the vicinity of the model. 
Subsequent correction comparisons will involve data from test 136, the early test of 
the six-inch-chord model and test 159, the test of the three-inch-chord model. 


M ref = 0.730, Re c = 10x10 6 
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Drag Rise Curves, C L = 0.3 


The following three figures show the Mach number correction in the form of 
uncorrected and corrected drag rise curves. The corrections are shown with and 
without the Murthy aspect ratio factor (ref. 6) on the Bam well- Se wall sidewall 
boundary layer term (refs. 3 and 5). The first figure below shows the comparison for 
Re c = 15 million at C L = 0.3. 


C L = 0.3, Re c = 15x10® 
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X10 


Drag Rise Curves, C L = 0.5 

Uncorrected and corrected drag rise curves are shown here for C L =0.5. 


C L ~ 0.5, Re c = 15x10® 
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OLX 


Drag Rise Curves, C L =0.7 


Uncorrected and corrected drag rise curves are shown here for C L =0.7. In all 
cases, the Bamwell-Sewall Mach correction is noticeably too large for the three-inch- 
chord test. The agreement is pretty good for all three cases with the Murthy aspect 
ratio factor included; this is taken as evidence that an aspect ratio factor should 
appear as part of a sidewall boundary layer approximation. 


C L = 0.7, Re c = 15x10® 
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Airfoil C p Plots 


The better correction to the Mach number due to the Murthy aspect ratio factor 
is also evident in comparisons of the free-air calculated pressure coefficient, shown by 
the vertical and diagonal crosses, and the experimental pressure coefficient 
renormalized with the corrected Mach number, shown by the asterisks. The figure 
shows the comparison for the three-inch-chord model at C L =0.37, M ref = 0.765, and 
Re c = 15 million. The shift in the Cp’s is eliminated by using the Murthy aspect ratio 
factor. Similar tendencies are found in the corrections for the six-inch-chord model 
but not to the same extent All subsequent corrections will be made with the Murthy 
aspect ratio factor included in the sidewall boundary layer approximation. 




Lift Curves and Error Parameter, M ref =0.60 

The following three figures show the results of applying the WIAC procedure 
to data for three Mach numbers and a Reynolds number of 15 million. The corrected 
and uncorrected lift curves are shown for two tests. In addition, e, the RMS matching 
error of the experimental and calculated airfoil surface velocity squared, is shown as 
an indication of the relative goodness’ of the corrections. As the error increases the 
corrections are deemed to be less trustworthy. The first figure shows the lift curve and 
the error for Mjg^O.60. 


M rS f = 0.600, Re c = 15x10® 



26 



Lilt Curves and Error Parameter, M ref =0.73 


This figure shows the lift curve and the error parameter for M^f-0.73, 
Re c = 15 million. It can be seen that the eiror parameter, e, becomes relatively much 
larger sooner with increasing a than was the case at ^^=0.60 shown on the previous 

page. 
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Lift Curves and Error Parameter, M ref = 0.765 

This figure shows the lift curve and the error parameter for 0.765, 
Re c = 15 million. It can be seen in these three figures that as the Mach number 
increases and the lift increases the error parameter also increases. This is due in part 
to the inability of the inviscid method to adequately model a flow condition greatly 
influenced by viscous and viscous/shock interaction phenomena. In addition, the 
present sidewall boundary layer/model pressure field interaction approximations may 
certainly become suspect at the higher transonic flow conditions. 
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CAST 10 Airfoil Data 


Comparisons of lift curve data for the French-built 18-cm (7.09 inch) chord 
CAST 10 model tested in both the NASA 0.3-m TCT and the ONERA/CERT T2 was 
recently given by Wolf and Ray (ref. 16). Both tunnels had adjusted wall test sections 
(AWTS) and both fixed and free transition results were given for Mr e f = 0.765 and 
Re c = 4 million. The curves shown in black by the squares andX’s on the figure below 
denote the fixed transition data. Lift curve data shown as open and closed circles on 
the figure are from a 6-inch-chord model tested in the 8- by 24-inch slotted wall test 
section (SWTS) of the NASA 0.3-m TCT with transition fixed at 7% chord. 
Uncorrected data are indicated by open symbols while the (second pass, 4-wall) 
WIAC data are given by the solid symbols. The GRUMFOIL free-air numerical 
results at the corrected conditions are denoted by an alternating dash-dot line when 
flow is attached (until very near the trailing edge) and a dotted line for separated 
fiow The value of c n appears to be larger for the slotted wall test section results. 

The corrected slotted wall data and GRUMFOIL results were taken from Gumbert 


and Newman (ref. 9). 

Mach = 0.765 ; Rc = 4 million 
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CAST 10 Airfoil Data 


Lift curve data shown by the solid lines is again that from the AWTS tunnels 
as described on the previous page. Data shown here by the open and closed circles 
is for “free” transition in the slotted wall test section of the TCT. Again the 
open symbols are uncorrected data, the filled symbols are WIAC data and the 
dashed curves are for GRUMFOIL free-air results. As pointed out in reference 9, 
the various GRUMFOIL results are for different transition locations (denoted n% 
at end of line) and it appears that the transition location in the tunnel tests 
is changing with lift level. The relative location of the curves in the present 
comparison indicates that the slotted wall test section appears to cause more 
premature transition than the adaptive wall test section. 



Angle of Attack, degrees 
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CAST 10 Airfoil Data 


Comparison of lift curve data for the Canadian-built 9-inch chord CAST 10 
model tested in the NASA 0.3-m TCT AWTS with that from the N AE 5-foot by 5-foot 
Blowdown Wind Tunnel with perforated top and bottom wall airfoil test section (15- 
bv 60-inch) were also given by Wolf and Ray (ref. 16). These results, shown by the 
X’s and the squares on the viewgraph below, are for transition fixed at 5 % chord at 
M , = 0 765 and Re c = 10 million. The Canadian data have been corrected for the 
top" and bottom perforated wall interference. Lift curve data shown as circles are 
from a 6-inch-chord model tested in the 8- by 24-inch slotted wall test section of the 
NAS A 0 3-m TCT with transition fixed (flagged symbols) and ’’free" (open symbols). 
The filled symbols represent the (second pass, 4- wall) WI AC data for free transition 
Broken line curves are again GRUMFOIL free-air results with transition denoted 
at the end of the curve. The shift in the angle of attack scale was simply due to 
different definitions for the zero angle of attack. 

Mach = 0.765 ; Rc = 10 million ; Transition Fixed 5% 



Angle of Attack, degrees wrt model chordline 



Angle of Attack, degrees wrt y=0 line 
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CAST 10 Airfoil Data 


Lift curve data shown by the X's and squares is again that from the NAE 
perforated and NASA AWTS tunnels as described on the previous page. Data shown 
here by the "plus" symbol arefor the 4- wall AWTS WIAC (ref. 4) applied to the NASA 
AWTS data by Mineck using the Murthy sidewall boundary layer option. The 
Navier-Stokes results denoted by the solid and dashed lines are due to Swanson et 
al and are discussed in the final talk of this workshop. At the higher lift levels for this 
Mach number, the Mach number corrections appear to be too large; apparently the 
subsonic wavy-wall solution invoked by Murthy (ref. 6) to approximately model the 
sidewall boundary layer effect is no longer valid for extensive supercritical flow and 
certainly not for large separated flow regions. This will also be discussed by 
Swanson. 


Mach 0.765 ; Rc - 10 million ; Transition Fixed- 5 % 
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CONCLUDING REMARKS 

PREMISE: All "airfoil tunnel" data contains 

some wall interference 


• Conclusions Concerning Data 

— Wall interference assessment must be made 

Wall interference corrections have been required to date 

— Corrections smaller for AWTS data than for SWTS data 
— SWBL influence can be significant at transonic 
high-lift conditions 

— Airfoil and tunnel-wall C p data required for TWINTN4 
— Transition location needs to be known 

• Conclusions Concerning WIAC 

— Transonic 4-wall approximations are required 
— Multiple passes needed to assess upstream flow angle 
— SWBL approx, needs to contain aspect ratio effect 
— Reasonable corrections seem to be obtained 
— Fairly easy to use 

— SWBL approx, may be inadequate for extensive 
supercritical flow 

— Interpretation of error parameter not yet established 
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